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ABSTRACT 

We study the statistical properties (such as shape and spin) of high- z halos likely 
hosting the hrst (PopIII) stars with cosmological simulations including detailed gas 
physics. In the redshift range considered (11 < z < 16) the average sphericity is 
(s) = 0.3 ±0.1, and for more than 90% of halos the triaxiality parameter is T < 
0.4, showing a clear preference for oblatcness over prolateness. Larger halos in the 
simulation tend to be both more spherical and prolate: we find s oc and T oc M^ T , 
with a s ~ 0.128 and = 0.276 at z = 11. The spin distributions of dark matter and 
gas are considerably different at z = 16, with the baryons rotating slower than the 
dark matter. At lower redshift, instead, the spin distributions of dark matter and gas 
track each other almost perfectly, as a consequence of a longer time interval available 
for momentum redistribution between the two components. The spin of both the gas 
and dark matter follows a lognormal distribution, with a mean value at z = 16 of 
(A) = 0.0184, virtually independent of halo mass. This is in good agreement with 
previous studies. Using the results of two feedback models (MT1 and MT2) by McKce 
& Tan (2008) and mapping our halo spin distribution into a PopIII IMF, we find 
that at high-z the IMF closely tracks the spin lognormal distribution. Depending on 
the feedback model, though, the distribution can be centered at w 65M© (MT1) or 
ps 140 Mq (MT2). At later times, model MT1 evolves into a bimodal distribution with 
a second prominent peak located at 35 — 40M Q as a result of the non-linear relation 
between rotation and halo mass. We conclude that the dark matter halo properties 
might be a key factor shaping the IMF of the first stars. 



1 INTRODUCTION 



The formation of first, metal-free (often referred to as 
PopIII) stars in the Universe represents a milestone dur- 
ing cosmic evolution, marking the e nd of the Dark Ages and 
producing the first heavy elements dCiardi fc Ferrara 20051; 
Yoshida et aHl200Sl : iBromm et al.l l2009; Br omm fc Yoshidal 



201ll ; lde Souza et al.ll201ll . l2012j ; [johnson et al1l2012l ). Thus, 
a key problem in physical cosmology is to understand the 
origin and evolution of such objects, born out of the pris- 
tine conditions leftover by the Big Bang. More specifically, 
the most urgent question concerns their Initial Mass Func- 
tion (IMF), which, despite its relevance, remains at best a 
poorly known quantity due to the lack of direct observations. 
Consequentl y, our knowledge is b ased mainly on theoretical 
models (e.g.. iDopcke et al]|2012t ). 

Until recently, studies based on the standard ACDM 
cosmological modej^ for structure formation predicted that 

Throughout the paper we assume a standard ACDM cosmologi- 
cal model, with current total-matter density parameter Q m = 0.3, 



the first stars formed when the age of the Universe was 
less than a few hundred million years, and that they were 
predominantly massive (lAbel et a! 1 120021 ; lOmukai fc Pallal 
l2003t lYoshida et alj|200d ) IClark et all l|201ll ). iGreif et al. 
1 20111 ') and IPrieto et aLl(|201lh have now performed cos- 
mological simulations using a sink particle technique to 
follow the evolution of a primordial protostellar accre- 
tion disk. They find that instead of forming a single 
massive object, the gas typically fragments into a num- 
ber of protostars with a range of different masses. How- 
ev er, high resolution rad iation-hydrodynamics simulations 
bv lHosokawa et all |201ll ) indicate a typical mass of PopIII 
stars of ~ 43Mp). Even mo re recently, similar results was 
found bv lStacv et ail (|2012i ). which suggests radiative feed- 
back will lower the final mass attainable by a PopIII star. 
With a estimative of 3OA/0 for a PopIII mass as a lower 



cosmological constant density parameter Q\ = 0.7, baryonic- 
mattcr density parameter Qf, = 0.04, expansion rate in units of 
100 km s _1 Mpc - , h = 0.7, spectral normalization erg = 0.9, 
and primordial spectral index n = 1. 
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limit. iGreif et al l (|2012h , using a different numerical scheme 
without need to insert sink particles, found that the kep- 
lerian disc around the primary protostar fragments into a 
number of secondary protostars, confirming previous results 
using sink particles. Clearly, these theoretical results are far 
from being conclusive, mostly due to the astonishing diffi- 
culty involved in simulating the large dynamical range re- 
quired and the complex physics involved. Note th at these 
studies, with few exceptions (e.g. iTurk et aI.ll2012T h largely 
neglect the possible effects of a magnetic field in the frag- 
mentation properties of the gas. 

In spite of this unsettled situation, a broad consensus 
exists on the fact that rotation of the protogalactic cloud 
is the key factor in determining the final outcome of the 
collapse. The importance of rotation has been fully appre- 
ciated after the fundamental paper bv lMcKee fc Tanl |200Sl. 
MT08), who studied the dependence of primordial protostar 
accretion in the presence of rotation (fostering the formation 
of an accretion disk) and radiative feedback from the pro- 
tostar. In their study MT08 concluded that the final stellar 
mass depends mainly on the entropy of the gas accreting 
from large radii, as well as its specific angular momentum. 
As the gas is bound in the gravitational potential of the dark 
matter halo, it follows that the angular momentum of the 
gas is linked to that of the parent halo. Hence, it is important 
to turn our attention to the properties of the dark matter 
halos that host the first stars to put the entire problem on 
a solid basis. 

The distribution of the angular momentum enters a 
broa d range of problems, from the halo m ass function itself 
(e.g. iMaccio et aT1l2007l : iDel Popoldl2009l), to the formation 
and evolution of central black h oles (e.g. lEisenstein fc Loebl 
1 19951 ; IVolonteri fc Reesl 120051 ; IVolonteril l20ld), to sem T 
analytical models of galaxy formation (e.g. lBensonll2012T ). 
Simulations of isolated dwarf galaxies suggest that angular 
momentum leads to more continuous star formation histo- 
ries than non-rotating cases by p reventing large-scale oscil- 
lations in the star formation rate ()Schroven et al.l[201ll ). For 
this reason, much effort has been done to explore the spin 
distribution of halos in different redshift and mass ranges, 
to characterize its probability distribution, as well as the 
dependence on mass, shape, merger rate and other halo 
properties (e.g.lBett et al.ll2007l ; iD'Onghia fc Navarrdl2007l ; 
lAntomiccio-Delogu et alj|20ld ). 

Another fundamental characteristic of halos is their 
shape. Spherical halos are very rare, and their collective 
propert ies cannot be appr oximated using spherical sym- 
metry (jAlkood et al.ll200d ). They are usually described in 
terms of ellipsoids characterized by three principal axis. 
Many authors have explored ways to estimate the shape of 
halos and the co rrelations with other halo properties (e.g. 
IZemp et al. Il201lh. Using a pure dark ma tter N-body simu- 
lations. Ijang-Condell fc Hernquistl (|200lh analyzed the pri- 
mordial halo characteristics and found no significant dif- 
ference between their results and simulations of large-scale 
structure formation at low redshift. iKazantzidis et all (|2004T l 
found that halos in cosmological simulations including gas 
cooling are considerably more spherical than those found in 
adiabatic simulations. This shows that the inclusion of de- 
tailed gas physics is fundamental, since the back-reaction 
effects of baryons on dark matter halos change their den- 
sity profiles as well as their mass distribution (|Cui et alj 



120121) . Observationally, weak gravitational lensing is proba- 
bly the m ost common approach to reconstruct the shape of 
halos (e.g. lCorless fc King||200Sl ; lBettll2012l ;l van Uitert et al.l 
l2012f ). The method does not depend on the presence of op- 
tical tracers and can be applied to a large range of scales. 

The purpose of the present paper is to explore the prop- 
erties (such as spin and shape) of the high-z halos likely 
to host the first stars, with an unprecedented inclusion of 
detailed gas physics. Although previous works that have 
inves tigated the correlation between halo parameters exist 
fe.g.. |jeeson-Daniel et al]|201ll ; [Skibba fc Maccicboill ). this 
is the first attempt to determine the characteristics and cor- 
relations of the low-mass end of the halo mass function in- 
cluding gas physics. 

As a final step, we embed our results into theoretical 
models of first star formation to unveil the links between 
dark matter halo properties and the PopIII IMF. 

The outline of this paper is as follows. In Sec. [2j 
we briefly describe the N-body/hydrodynamical simulations 
used to derive the halo properties; Sec. [3] describes the 
methodology used to calculate the quantities of interest 
(halo spin and shape). In Sec. 2] we show the results and 
provide useful fits for correlations between different quan- 
tities. Finally, Sec. [5] discusses the possible implications for 
the PopIII IMF. Sec. [5] contains a summary of the results. 



2 SIMULATIONS 

We analyze the outpu t of the N-b o dy/hy drodynamical sim 



ulations described in iMaio et~aT 



>m L /nyc 
(l20ld). 



which were per- 



formed using the Gadget-2 code 
lations include the evolution of e~ 
He ++ , H 2 , Hj" Di D+, HD, HeH 



Springelll2005l ). 

H H+, H~, He He+ 
(jYoshida et al.l 120031 ; 
IMaio et al ll200d . 120071 , l2009h . PopIII and PopI I /I sta r for- 
mation and metal pollution (|Tornatore et al.l 20071) . gas 
cooling from resonant and fine-struc ture lines (Maio et al. 
20071 . I2OO9I ) and feedback effects ([Springel fc Hernouist 
20031 ). There are evidences for the existence of a criti- 
cal metallicity, Z cr it, which allow the tran sition between 
PopII/I and Pop H I star formation modes (|Omukail |2000| : 
iBromm et al.ll200ll ). The transition from the PopIII to the 
PopII/I regime is determined by the value of the gas metal- 
licity, Z, compared to the critical value Z cr u = W~ 4 Zq. The 
value of this minimum level is very uncertain, but is likely to 
be between 10~ 6 and 10 _4 Mq and can be strongly depen- 
dent of the efficie ncy of dust formation in fi rst-generation 
supernova ejecta (|Schneider et al.l [2003, 200^) and the fine- 
structure line c ooling of singly-ionized carbon or neutral 
atomic oxygen (|Bromm fc Loebl |2003l ). If Z < Z cr i t a 
Salpeter IMF (i.e. with a slope of -1.35) is assumed in the 
mass range 100-500 M© , otherwise the same Salpeter slope is 
adopted in the mass range 0.1-100 Mq. The chemical model 
follows the detailed stellar evolution of each SPH particle. At 
every time-step, the abundances of various heavy elements 
(C, O, Mg, S, Si and Fe) are consistently derived, according 
to the lifetimes and metal yields of the stars in the given 
mass range. 

The cosmological field is sampled at redshift z = 100, 
with dark-matter and baryonic-matter species. We consider 
snapshots in the range 11 < z < 16, within a cubic vol- 
ume of comoving side 1 Mpc h~ x and 320 3 particles per gas 
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and dark matter species, corresponding to particle masses 
of 116 Mq/i -1 and 755 MqIi -1 , respectively. The identifica- 
tion of the simulated objects is done by applying a Friends 
of Friends (FoF) technique; substructures are identified by 
using a SubFind algorithm l|Dolag et al.l l2009h . which dis- 
criminates among bound and non-bound particles. For more 
details on the simulation s we refer the reader to the original 
paper (|Maio et alj|2010h . 



3 HALO PROPERTIES DERIVATION 

In the following we describe the method used to derive the 
halo properties of interest here, as its shape and spin, along 
with a number of ancillary quantities defined below. 



a weight factor, w(r), can be introduced 

S ^-N^ d 2 ■ (5) 

%= i 1 

A critical review o f the different approaches can be found in 
IZemp et al.l (|201ll ). They explore six methods, differing by 
the integration volume and the choice of the weight func- 
tions, w(r) — l,w(r) — l/r 2 ,w(r) = l/df. They conclude 
that using weights can introduce a systematic bias in the 
measured axis ratios, in addition to blurring the physical 
interpretation of the shape tensor. Thus, using w(r) = 1 
and integrating over the enclosed ellipsoidal volume appears 
to be the most unbiased method. This choice is also pre- 
ferred for halos with lower number of particles and if the 
main interest is in deriving the global rather than local (i.e. 
as a function of distance) shape. Hereafter, all calculations 
are done using eq. [2] integrated over the enclosed ellipsoidal 
volume. 



3.1 Shape 

The halo shape is estimated based on its mass distribution, 
which can be direct ly derived using the eigenvalues of the in- 
ertia tensor X (e .g., Ispringel et alj|2004 lAllgood et alj|2006l ; 
iBett et alj|2007h . 



Xj k = y]mi(rj5j,k 



'i,jii,k ) 



(1) 



where rj and m; are the position vector and mass of the i-th 
particle, 5j t k is the Kronecker delta and the sum is performed 
over the total number of particles inside the halo, N. 

An alternative way to measure the shape is by using 
the second moment of the mass distribution, i.e. the shape 
tensor 



1 N 
S » ~E' 



(2) 



As the eigenvalues of S and X are the same, the two defini- 
tions are totally equivalent for halo shape studies. We then 
restrict our discussion to the shape tensor. 

The eigenvalues of the diagonalised shape tensor define 
an ellipsoid, which represents the equivalent homogeneous 
shape of the halo in terms of the principal axis ratios, with 
the convention a ^ b ^ c. It is customary to refer to the axis 
ratios in terms of sphericity, s — c/a (with s = meaning 
aspherical and s = 1 spherical), oblateness, q = b/a and 
prolateness, p = c/b. With these definitions, the triaxiality 
parameter can be conveniently written as 



T 



1 



(3) 



hence a prolate (oblate) halo has T = 1 (T = 0). 

A slightly di fferent way to calcula te the shape tensor 
was introduced by lAllgood et all (|2006l ). 



1 N 

-Y 



N 



(4) 



where d 2 = x 2 + y 2 /q 2 + zf / s 2 is the elliptical distance in 
the eigenvector coordinate system from the center to the i- 
th particle and q and s are found iteratively. More generally, 



3.2 Spin 

The spin parameter is a measure of the amount of coherent 
rotation in a system compared to random motions. For a 
spherical object, it corresponds approximately to the ratio 
of its own angular velocity to the angular velocity needed for 
it to be supported against gravity solely by rotation (see e.g. 
Padmanabhan 1993). The halo spin can be characterized by 
a dimensionless parameter A, 



A : 



J\E\V 2 



(6) 



GM 5 /2 ' 

where J, E, and M are the total angular momentum, energy, 
and mass of the system, and G is the gravitational constant. 
The (specific) angular momentum per unit mass is 



1 

vE 1 



N 



with n and Vi being the position and velocity of i-th particle 
relative to the halo center and halo center of momentum, 
respectively; is again the total number of particles inside 
the halo. The kinetic, Ek, and potential, E p , energies are 
calculated on-the-fly during the simulation as 



E k 



1 N 



E v — 



f N 2 




(- Gml ) 


\N 2 - 


- Nu) 





JV«-1 N u 

E E 

»=1 3=»+l 



where r\ is the softening length and W(u) is the softening 
kernel. If the halo contains more than 1000 particles, the po- 
tential is calculated using N u — 1000 randomly selected pa r- 
ticles (for more details see lSpringeHl2005l ; IBett et al.ll2007l ). 



4 RESULTS 

4.1 Shape distribution 

The shape dependence on halo mass has been considered 
previously by several authors. Despite they all agree on the 
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Figure 1. Sphericity s (upper panels), and triaxiality T (lower 
panels), as a function of total halo mass for halos at z = 16 
(upper figure) and z = 11 (lower). The contour levels represent 
50%, 75%, 90%, 95% and 97.5% of the sample. 



non-spherical nature of halos, overall conclusions can be dif- 
ferent depending on the assumptions made to define ha- 
los, the methods to measure shapes or t h e incl usion of gas 
physics (|Allgood et alj|2006t ). iBett et all (|2007l ) found that 
mor e massive halos tend to b e less spherical and more pro- 
late; iKazantzidis et all f2004) noticed that halos formed in 
simulations with gas cooling are more spherical than halos 
in adiabatic simulations. 

To estimate the halo shape, we only use halos with more 
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Figure 2. Spin distribution of the total (dark matter + gas) halo 
(black dashed lines), the dark matter (grey dot-dashed lines) and 
the gas (solid lines). The upper and lower panels refer to z = 16 
and 11, respectively. 



than 100 particles (gas + dark matter), equivalent to a total 
mass of ~ 10 4-5 Mq. The shape of the halos is described in 
terms of sphericity, s, and triaxiality, T, as defined in Sec. 
3.11 The probability contour levels of s and T as a function 
of the total halo mass, Mh, are shown in Fig. [T]for z = 16 
(upper panel) and z = 11 (lower panel). 

In the entire redshift range considered (11 < z < 16) 
the average sphericity is (s) = 0.3 ± 0.1, and for more than 
90% of halos T < 0.4, showing a clear preference for oblate- 
ness over prolateness. This is markedly different from 2 = 
halos that tend to be more prolate a nd spherical: for ex - 
ample, pure collisionless simulations (|Allgood et al.l I2006T ) 
found (s) ~ 0.6±0.1 for galaxy mass halos ~ 10 12 /i _1 M Q . 
Note that constraints from observations of the Sagittarius 
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tidal streams give best-fit parameters s ~ 0.67, q ~ 0.83, 
and T ~ 0.56, in agreement with the Galactic model by 
lLaw et all (|2009h . 

Using a non-linear last square method computed with 
the Gauss-Newton algorithm, and approximating s oc M^ s 
and T oc M^ T , we find a a « 0.147 and a T = 0.285 at z = 
16, and a a « 0.128 and qt = 0.276 at z = 11. The fits are 
given belowj, 



<«(*)> = C 



/i-'M, 



(9) 



-4.1082 + 0.9834z - 0.07542 2 + 0.001912 3 , 



and 

O.T 

Ct 



= 20.143 - 4.9772 



0.3842 2 - 0.00979z 3 , 



= c 



h 



h~ 1 M Q , 
-9.7534 + 2.29242 



(10) 



0.17272 2 +0.004292 3 , 



= 51.720 - 12.5002 + 0.9462^ 



0.02362 J 



The most massive halos in the simulation (M/j = 
1O 6 - 5_7 M ) tend to be more spherical and prolate than 
smaller ones, with weak variation of the mass dependence 
with redshift. 
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Figure 3. Cumulative distribution function of dark matter spin 
(thick blue solid line) compared with theoretical distributions at 
z = 16. The theoretical distributions arc: lognormal (black solid 
line), beta (red dashed), Weibull (orange dotted) and Gamma 
(purple dot-dashed). 



4.2 Spin distribution 

The spin distribution of the total (dark matter + gas) mass, 
dark matter mass and gas mass is shown in Fig.[2]for z = 16 
(top panel) and z = 11 (low panel). The curve is smoothed 
using a kernel density estimator for a sample of n elements, 



f(x,h s ) = 



nh a {x 



J2 K 



h s {x) 



(11) 



with a Gaussian kernel K and an adaptive bandwidth h s - 
The distributions of dark matter and gas are considerably 
different at high redshift (z = 16), with the baryons rotat- 
ing slower than the dark matter, which gives the dominant 
contribution to the total spin. At lower redshift, instead, 
the spin distributions of dark matter and gas track each 
other almost perfectly, as a consequence of a longer time 
interval available for momentum redistribution between the 
two components. It is important to notice that the com- 
parisons are done by collecting different dark matter and 
gas particles of the same halo, thus fully accounting for the 
back-reaction of baryons on the parent dark matter dis- 
tribution. It is very common to fit the halo spin distribu- 
tion with a lognormal f unction jBailin fc Steinmetzl 120051 : 
iDavis fc Nataiaiarj20ld ) . However lBett et all (|2007l l found 



deviations from such functional form when studying a large 
(> 10 6 ) number of halos in the Millennium simulation. Both 
the lognormal and Weibull models can be used quite ef- 
fectively to analyze skewed data sets. Although these two 
models may provide similar data fit for moderate sample 
sizes, the inferences based on the model will often involve 
tail probabilities, where the effect of the model assumptions 
is very critical. This makes it important to use a quantitative 



2 Note that all fits hereafter are valid only within the redshift 
range of the simulation, unless explicitly stated. 



diagnostic to quantify the best distribution to use. To do so, 
we test four classical distributions: lognormal, Gamma, beta 
and Weibull, whose shapes are given by: 

1 



f L (x;n,cr) = 

fr(x-k,6) = 

fp{x;a,p) = 

fw(x;k,\) = 



-(lnx-(i) /2(T 2 



xa 
1 



2tt 
1 



-1 (1- 



x)/e 



!> + /?) 



r(a)r(/3)x° 

k /in*- 1 
A (a) ' 



-Hi- 



(12) 
(13) 
(14) 
(15) 



In order to choose the most suitable distribution, we fit the 
spin distribution of halos in our redshift range. Then, we use 
a Maximum Likelihood test to justify our choice. We obtain 
the following redshift-averaged values of the reduced chi- 
square: (xl,x},xl,Xw) = (2.89±0.94, 11.95±4.85, 11.38± 
4.85, 13.53 ±4.34). In Fig. 03 we plot the cumulative distri- 
bution function for the four distributions above and the one 
obtained from the simulations. It is clear that the lognormal 
fits the data at best. 

The distribution of spin parameter can be written as 



P(A) 



1 



A<ToV27T 



exp 



In 2 (A/Ao) 



2ol 



(16) 



where Ao is the location parameter and 00 is the shape pa- 
rameter of the distribution. The best fit values for mean and 
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variance, a 2 , as a function of redshift arfjfl 
(X(z)} = 1.315-3.681 x 10~ 1 z + 3.903 x 1CT 



1.831 x l(T 3 .z 3 +3.206 x 10~V, 



and 



a(z) = 0.1754 - 5.246 x 10 z + 5.871 x 10" 
- 2.909 x 10" 4 2 3 + 5.385 x 10~ 6 z 4 . 



(17) 



(18) 



While it is well known that the dependence of spin 
on halo masses is relatively weak ()Macci6 et all 120071 : 
iBett et al]|2007l '). it has not been verified yet whether this 
holds also for small objects. Therefore, analogously to the 
shape distribution, in Fig. 0] we show the contour levels of 
the halo spin as a function of total halo mass at z = 16 and 
11. 

Assuming a power law A oc M£ , best fit becomes 



<A(*)> = Ca 



M h 



a x = -1.64 x 10 2 + 4.97 x lO 1 ^ 
+ 2.78 x 10" 



(19) 



5.60z 



-1 3 
z 



5.16 x 10" 



3 4 
z , 



Ca = 4.056 - 0.993^ + 0.045z 



0.004z J 



The slope ct\ evolves from —0.023 at z — 16, to 0.012 at 
z = 11; both values are consistent with 0, indicating that 
the spin parameter distribution is essentially independent 
on halo mass also at the very high redshifts considered here. 



5 IMPLICATIONS FOR THE POPIII IMF 

In the following we discuss the possibility to connect the host 
halo properties to the typical mass scale of the collapsed pro- 
tostar. It is important to keep in mind that all results pre- 
sented here rely on the best available semi-analytical model 
to translate the angular momentum distribution into a corre- 
sponding stellar IMF. Despite such limitation, our approach 
provides the best way to statistically analyze a large sample 
of simulated halos, which could not be performed otherwise. 
For a detailed study of the properties of the rotation and 
st ructure of PopIII st ars we refers the r eader for the works 
of iGreif et all (|2012T ); IStacv et al.l (|2012T ). They resolve four 
minihalos down to scales as small as 0.05 Rq. They found 
little evidence of correlation between the properties of each 
host minihalo and the spin of its largest protostar or the 
total number of protostars formed in the minihalo. However 
due the low number of halos probed, it's difficult reach a 
statistical significant conclusion. 

The IMF is usually defined as a segmented power-law or 
a log- normal type mass distribution (|Kroupall200"ll ; IChabrierl 
120031 ). Both definitions are correlated by 

dN oc m* a drn- k or dN oc m* d(log m*), (20) 

where T = -(a - 1) |Bonnell et al.l l2007f ) and N is the 
number of stars with masses in the range m, to m* + dm*. 
The Salpeter slope jSalpeterl ll955T l is given by a = 2.35, 
or r = —1.35. As discussed in the Introduction, the PopIII 
IMF is actually unknown, but there are hints that it could 
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Figure 4. Spin distribution as a function of the total mass for 
halos at z = 16 (top figure) and 11 (bottom). The contour levels 
in the upper panels represent 50%, 75%, 90%, 95% and 97.5% of 
the sample, while in the lower panels the median values per bin of 
mass (log Mfj = 0.5) is shown. The width of boxes is proportional 
to the square root of the number of halos within each bin and the 
whiskers extend to the most extreme data point, which is within 
the 50% interquartile range, i.e. the difference between the largest 
and smallest values in the middle 50% of the dataset. 



3 The mean and variance of lognormal distribution are given by 
(A) = e A o+ CT o/ 2 and a 2 = e CT o+ 2A o( e ^ _ iy 
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have been biased towards more massive stars compared to 
that of PopII/I stars. 

One of the important ingredients for the determination 
of the IMF is the typical rotation of the gas. Extreme ultravi- 
olet radiation from the protostar can ionize infalling neutral 
gas, creating an H II region whose expansion reduces signifi- 
cantly the accretion of gas onto the star. The radiation would 
also destroy molecules and inhibit star formation in the 
surroundings, by affecti ng accretion onto more distant pro- 
toste l lar cores (see e.g. |Ricotti et al.ll2002l; lAhn fc Shapirol 



120071 : IWhalen fc Normanll200Sl ; IPetkova fc Maioll2012l ). The 

expansion is facilitated by protostellar cores with higher ro- 
tation, as the density in the polar directions tends to be 
lower; hence the final stellar mass is expected to increase 
in slowly rotating cores where the effects of feedback are 
quenched. 

In the following we will make a simple estimate for the 
expected PopIII IMF, based on a combination of results 
from our simulations and the semi-analytic prescription de- 
scribed in MT08. In addition to studying the dependence 
of the final stellar mass on several radiative feedback pro- 
cesses such as photodissociation of H2, Lyman-a radiation 
pressure, formation and expansion of H II regions, and disk 
photo-evaporation, the authors investigate the role of gas ro- 
tation, which modifies the density distribution in the vicinity 
of the star. In their fiducial model they assume that accre- 
tion is halted when the disk photo-evaporation rate exceeds 
the accretion rate onto the star-disk system; at that stage 
the typical PopIII stellar mass is ~ 137M0. 

Such value depends strongly on the assumed rotation 
rate of the protostellar core, which in their case was taken 
to be half of the keplerian velocity as sugges ted by AMR sim- 
ulations focusing on a single parent halo (lAbel et al.l |2002| ; 
iBromm et al.ll2002l ; lo T Shea fc Normanll2007f ). However, high 
rotational speeds might not be the rule, as we h ave seen from 
our previous analysis. IHosokawa et al.l (|201lh have found 
qualitatively similar results by means of a 2D hydrodynamic, 
radiative transfer simulation, showing that the mass accre- 
tion is shut down due to the dynamical expansion of the 
Hn region and the photo-evaporation of the circumstellar 
disk. However, their simulations found systematically lower 
final masses than MT0 8. The likely explanation suggested by 
IHosokawa et al.l (|201ll ) for this discrepancy comes from the 
different model of the stellar feedback. In addition, MT08 
assume that after the formation of the Hn region mass ac- 
cretion onto the disk still continues from regions shaded by 
the disk, which is not in agreeme nt with the simulation pic- 
ture from IHosokawa et al.l (|201ll ). 

In the following, we consider the influence of rotation 
in two different scenarios developed by MT08. In the first 
(hereafter MT1), the accretion is reduced by the expansion 
of the H 11 region around the prot ostar and sets a typica l 
mass similar to the values found by IHosokawa et all (|201lf ). 
In the second (hereafter MT2), the accretion efficiency is 
reduced as the H 11 region expands, however, accretion is 
allowed to continue from directions that are shadowed by 
the disk photosphere. The accretion stops when the photo- 
evaporation rate exceeds the accretion rate onto the star- 
disk system and sets a larger typical mass scale. The com- 
parison between our results and MT08 can be easily done 
due to the one-to-one relationship between fk ep defined be- 
low and the stellar mass scale, assuming other parameters 



fixed (see Fig. 5, 10 and Table 1 in MT08, where the values 
to interpolate the one-to-one relationship were taken). 

The angular momentum of the gas accreting onto the 
star-disk system can be characterized by its keplerian pa- 
rameter: 



f kep (M n ) = V^/V l 



kep,i ■> 



(21) 



where V^,j and Vk ep ,i are the rotational and keplerian ve- 
locity, respectively, as a function of the total mass enclosed 
within a radius r;. Here, is the position vector of the z-th 
particle relative to the halo centre. The characteristic veloc- 
ities above are defined as: 



V^,i = ji/n, Vke P ,i = \/GM(ri)/ri, 



(22) 



where the specific angular momentum, ji, of the i-th particle 
is averaged over the spherical shell whose radius is Vi. 

The value of fke P averaged over all particles within a 
halo can be expressed as 



(fkep) 



J_ -sp (ji/n) 

N ^ 



kep.i 



(23) 



According to iMcKee fc Tanl (|2008l) . for f kep > 0.25 little 
difference is observed in the final stellar mass, which is set 
by the balance between the (inner) disk-shadowed accretion 
and mass loss due to photo-evaporation. For smaller rota- 
tion parameters (fke P < 0.125) instead, the mass scale at 
which accretion is halted strongly depends on the H II re- 
gion breakout. 

To calculate the (fkep) distribution, we use the same 
procedure described in Sec. 14.21 and we show it in Fig. [5] at 
2 = 16 and 11. The overall shape of the distribution remains 
qualitatively similar at different epochs, however the mean 
value increases towards lower redshifts, from 0.26 at z = 16 
to 0.61 at z = 11. 

The typical rotational velocity is in general below the 
required velocity for rotational suppo rt, in agreement with 
previous calculations ()Abel et al.ll2002l ). Using the above dis- 
tributions we can translate the rotational velocity into a 
PopIII typical mass using models MT1 and MT2. The re- 
sult of this exercise is shown in Fig. [6] there are several 
interesting features that we can deduce from here. 

First, at high redshift the IMF tends to closely track the 
lognormal distribution imprinted by the rotation properties 
of the halos. Depending on the feedback model, though, the 
distribution can be centered at w 65M (MT1) or w 140M Q 
(MT2). At later times, model MT1 tends to evolve into a 
bimodal distribution with a second prominent peak located 
at 35 — 40M© in addition to the initial one. The bimodality 
comes from the non-linear connection between rotation and 
mass scale. For values of fk ep > 0.25 the rotation has a weak 
influence on the final stellar mass. As the redshift decreases, 
so does the width around the second peak of the stellar mass 
distribution, because the majority of halos have higher val- 
ues of /feep- A peak at m* ~ 65Mq is still present due to the 
slow rotation tail of the (fke P ) distribution. Model MT2 in- 
stead shows a much more gradual and moderate shift of the 
peak towards lower masses, accompanied by an increasingly 
narrower width. 

Thus, it seems that the (fk ep ) distribution and shift 
with redshift, governed by the angular momentum evolution 
of the halos, has an extremely strong influence on the PopIII 



8 R. S. de Souza, B. Ciardi, U. Maio, A. Ferrara 




Figure 5. Distribution of log(/fc ep ) at z 
11 (dashed green). 



=16 (solid blue line) and 



q 
d 



IMF, at least as long as we assume the MT08 feedback mod- 
els to be correct. However, it is important to emphasize that 
the gas rotation is also affected by thermal heating from 
feedback mechanisms, which are expected to be stronger for 
PopIII stars and in low mass halos, a s the ones we are deal- 
ing with. In addi t ion, s imulations fe.g. lTornatore et alj2007t 
iMaio et al.ll2010l .l 201 ll ) show that it is common to have mul- 
tiple star formation sites within the same halo, with a com- 
bination of PopIII and PopII/I stars. This means that it is 
not straightforward nor trivial to assign a single PopIII IMF 
to a halo. 

With the above caveats in mind, we come to the some- 
what surprising conclusion that, although on a protostellar 
basis radiative feedback acting on baryons might be the key 
factor in determining the mass of the first stars, it is the an- 
gular momentum distribution of the dark m atter halos that 
contr ols the build-up of the IMF (see also ISchroven et al.l 
2011). This process might work in the simple way outlined 
here as long as there is a one-to-one correlation between 
the halo and the protostellar core angular momentum and 
it might break down in larger galaxies in which momentum 
is dissipated via tidal torques and/or shocks arising from 
the interaction among different cores or galactic-scale dy- 
namical instabilities. We reiterate that all our calculations 
include the back-reaction of baryons on the dark matter, 
and hence they should provide a robust description of the 
total matter dynamics in a halo. 



6 SUMMARY 

Our study, following the evolution of dark matter and bary- 
onic physics in cosmological simulations, makes possible to 
study the statistical properties of the high-z, low mass halos 
that likely hosted the first stars. In addition such simula- 
tions include a large number of physical processes (PopIII 
and PopII/I star formation, metal enrichment, gas cooling 
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Figure 6. PopIII typical mass distribution for model MT1 (top 
panel) and MT2 (bottom panel) at redshift z =16 (blue solid 
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12 (brown long-dashed), and 11 (green long dashed-dot). 



from resonant and fine-structure lines and feedback effects) 
and a detailed chemical network following the abundances 
of key species (e~, H, H+, HT, He, He+, He ++ , H 2 , Hj, D, 
D+, HD, HeH+). 

In this work we have mostly concentrated on the sta- 
tistical analysis of two important halo properties, i.e. spin 
and shape. As these parameters, governing the overall evo- 
lution of protostellar cloud collapse, are predicted by mod- 
ern PopIII formation theories to be related to the mass of 
the first stars forming in these systems, we then discuss the 
implications of our findings for the PopIII IMF. Our main 
results can be summarized as follows: 



In the entire redshift range considered (11 < z < 16) 
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the average sphericity is (s) = 0.3 ± 0.1, and for more than 
90% of halos T < 0.4, showing a clear preference for oblate- 
ness over prolateness, contrary to what found at z = 0. 

• Larger halos in the simulation tend to be both more 
spherical and prolate: we find s oc M^ s and T oc M^ T , with 
a a « 0.128 and a T = 0.276 at z = 11. 

• The spin distributions of dark matter and gas are con- 
siderably different at z = 16, with the baryons rotating 
slower than the dark matter (giving the dominant contri- 
bution to the total spin). At lower redshift, instead, the spin 
distributions of dark matter and gas track each other almost 
perfectly, as a consequence of a longer time interval available 
for momentum redistribution between the two components. 

• The spin distribution for both gas and dark matter in- 
side the simulated small halos can be well represented by 
a lognormal function, with mean and variance at z = 16 
of 0.0184 and 0.000391, virtually independent on halo mass 
and in good qualitative agreement with previous results. The 
mean value of spin parameters is also in agreement within 1- 
a with the median value found bv lJang-Condell fc Hernquistl 
i|200lf ) for a lognormal distribution in their study of small- 
scale structures at high-z using a pure dark matter N-body 
simulation. 

• According to most recent theories of PopIII star for- 
mation, rotation is the key factor in determining their final 
mass. Using the results of two feedback models (MT1 and 
MT2) by McKee & Tan (2008) and mapping our halo spin 
distribution into a PopIII IMF, we find that at high z the 
IMF tends to closely track the spin lognormal distribution; 
depending on the feedback model, though, the distribution 
can be centered at « 65M Q (MT1) or « 140M Q (MT2). At 
later times, model MT1 tends to evolve into a bimodal distri- 
bution with a second prominent peak located at 35 — 40Mq , 
as a result of the non-linear relation between rotation and 
halo mass. 



Po pIII IMF is still highly d e bated 



While the 

dStacv et al l l20ld; Idark et al.1 l201ll: [Greif et al 
iPrieto et al.ll201ll ; Iflosokawa et all 1201 ll ; iGreif et al 



the present study offers an intriguing indication that the 
IMF of the first stars might be tied and controlled by the 
properties of their parent halos, thus linking in a novel 
way large scale structure and early star formation. If this 
is indeed the case, our suggestion could lead to clear and 
testable predictions (e.g. PISN rates, abundance of pure 
PopIII galaxies, metal abundance patterns in the IGM and 
low-mass stars to mention a few) for the number, properties 
and cosmic evolution of these pristine stellar systems. 
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